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Abstract 

We explore similarities in the luminosity distribution of early type galaxies and the mass profiles of ACDM 
halos. The spatial structure of these systems may be accurately described by a simple law where the logarithmic 
slope of the projected density is a power law of radius; the Sersic law. We show that this law provides a 
significantly better fit than a three-parameter generalization of the NFW profile and derive the best-fitting 
Sersic parameters for a set of high-resolution ACDM halos spanning a wide range in mass. The mean Sersic 
n values are 3.0 for dwarf- and galaxy-sized halos and 2.4 for cluster-sized halos, similar to the values that 
characterize luminous elliptical galaxies. We discuss possible reasons why the same law should describe dark 
and luminous systems that span a range of over seven decades in mass. 
Subject headings: 



1. INTRODUCTION 
The lS6r^lfT96Ri) law. 

ln(E/E,) = -/7(X 1 /"-i; 



(1) 



relating the 2D (projected or surface) density, E, and the di- 
mensionless radius, X = R/R e , is often fit to the luminosity 
profiles of elliptical galaxies and to the bulges of disk galax- 
ies. The parameters of the fit include the Sersic index, n, 
as well as the constant, b, which is normally chosen so that 
R e is the radius containi ng one-half of the projected light; 
b = b(n) « 2n - 0.324 ( Ciotti & Berti nll999l) . 

In a recent series of papers, A. Graham and co-workers 
have shown that the Sersic law provides a remarkably good 
fit to the luminosity profiles of stellar spheroid s, from dE 
galaxi es to the most luminous ellipticals ( Graham 2001, 2002; 
Graha m & Guzmanll2003t IGraham et alJl2003t iTrujillo et all 
2004). The fits apply over 2-3 decades in radius, and often 
extend down to the innermost resolvable radius. Deviations 
from the best-fitting Sersic law are typically of order 0.05 
magnitudes rms. The Sersic index n is found to correlate well 
with galaxy absolute magnitude, 

log 10 fl«-0.106M B -1.52 (2) 

JGraham & Guzmanl 120031) . and also with other structural 
parameters like R P and E f JCaon. Capaccioli & D'Onofriol 
1 19931 IGraham & Guzmal] 120031) . Setting n = 4 gives the 
Ide Vaucouleursl i ll 9481) law, which is a good fit to luminous 
elliptical galaxies, and n = 1 is the exponential law, which 
reproduces well the luminosity profiles of dwarf ellipticals. 

There are some known limitations to the applicability of 
equation Q to the very central regions of some galax- 
ies. In particular, Sersic's law fails to represent adequately 
the very central profiles of elliptical galaxies with cores; 
the pointlike nuclei of some dE galaxies; and the steep 
power-law density cusps observed in the inner few parsecs 
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of nearby galaxies like M32 and the bulge of the Milky 
Way. The origin of these features is not well understood, 
but it is likely that they are the result of dynamical pro- 
cesses, possibly involving single or multiple black holes, 
which act to modif y the pre-existing Sersic profile in the 
inner most regions {v an der Marel 1999; Milosavlievic et all 
1 20021 iRavi ndranath. Ho & Filippenko 2002; Graham l2004F 
iMerritt et alJ2004HPreto. Merritt & SpurzerrJ2004l) . 

The density profiles of the dark matter halos formed in 
A/-body simulations of hierarchical clustering have tradition- 
ally been fit to a rather dif ferent class of functions, es- 
sentially broken power l aws jNavarro. Frenk & Whitdfl996l 
ll997t Mo^r^^lJll^99l) . However, the m ost recent simula- 
tions JPower et alJl2003b IReed et alJl2004l) suggest that halo 
density profiles are better re presented by a function with a 
continuously-varying slope. Navarro et al. (2004) proposed 
the fitting function 



dlnp/dlnr = -2 (r/r_ 2 ) a 



(3) 



where r_2 is the radius at which the logarithmic slope of the 
space density is —2, and a is a parameter describing the de- 
gree of variation of the slope. The corresponding density pro- 
file is 



ln(p/P-2) = -(2/oc)(* a -l) 



(4) 



with x = r/r-2- Remarkably, this is precisely the same func- 
tional form as equation - with the difference that Navarro 
et al. fit equation l|4} to the space density of dark matter ha- 
los, while equation Q applies to the projected densities of 
galaxies. 

Nevertheless the connection is intriguing and a number of 
questions spring to mind. Does the Sersic profile fit the sur- 
face density profiles of dark matter halos as well as it fits 
galaxies? We will show here (§3) that the answer is "yes": the 
same fitting function provides an equally good description of 
the projected densities of both dark and luminous spheroids. 
In §4 we ask whether it is most appropriate to fit the Sersic 
law to the space or projected densities of dark matter halos, 
and whether these functions are better fits than other three- 
parameter functions. §5 contains some speculations about 
why a single density law should describe dark and luminous 
systems over such a wide range in mass. 

2. METHOD 
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FIG. 1 . — (a) Nonparametric estimates of the surface density profiles of the 
19 halo models. Profiles of the D (C) models have been shifted downward 
(upward) by 0.75 in the logarithm, (b) Deviations of the best-fitting Sersic 
model from £(R). Fitting parameters are given in Table 1. 



We constructed nonparametric estimates of the space and 
projected den s ity pr ofiles of the 19 ACDM halo models in 
iNavarro et alJ ( 120041) . and compared them with a number of 
fitting functions, including the Sersic law, equation Q; the 
deprojected Sersic law £ ( /(r), defined as the spherical den- 
sity law whose spatial projection is E(T?); and a generalized, 
three-parameter NFW (1996, 1997) profile, which may be ex- 
pressed as 



dlnp/dlnx = -(y+3x)/(l + x) 



(5) 



with x = r/r s . The NFW profile has y = 1 and r s = r_2; the 
iMoore et al.l (TT999) profile has a similar functional form with 
inner slope y= 1.5. 

Details of the numerical simulations are given in 
INavarro eta D <|2004). Four halos are "dwarf" sized (M w 
10 10 M Q ), seven are "galaxy" sized (M « 1O 12 M ), and eight 
are "cluster" sized (M as 10 15 M Q ). We adopt the notation of 
that paper (D= dwarf, G=galaxy, C=cluster) in what follows. 

Nonparametric estimates of the space and projected den- 
sity profiles, p(r) and £(/?), were constru cted using the s pher- 
ically symmetriz ed kernels defined by iMerritt & Tre mblav 
( 1994) (see e.g. iReed et alJ (120041) . Appendix A). Each TV- 
body point was replaced by a kernel of the form 



K p (r,ri,hi) ■■ 



1 



2(27t) 3 /2 \hj 



-(^)/2^ sinh ( rr ./,^ 



TABLE 1 

Model fits to the halo density profiles. 
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K z (R,R„hi) = -L e -( R - +R2 )l {RRi/hl) 
2k 



(7) 



with hi the width of the kernel associated with the /th parti- 
cle and To the modified Bessel function. The projected radii 
Rj were obtained from the TV-body radii r, by assigning each 
particle a random position on the sphere of radius r,-. Density 
estimates were computed on a grid of 100 radial points spaced 
logarithmically from r conv to r2oo (thes e radii are defined be- 
low). We followed standard practice (Silverman 119861) and 
first computed a pilot estimate of the density via a nearest- 
neighbor scheme, then allowed the /i, to vary as a power 8 of 
this pilot density. 

When fitting one of the parametric functions defined above 
to p or £, we computed the density estimates on a grid in 
(ho, 8) (Iiq is the geometric mean of the hj) to see which choice 
of kernel parameters minimized the residual of the fit; typ- 
ically there was a broad range of (/io,8) values over which 
the best-fit parameters and their residuals were nearly con- 
stant. The residual was defined as the rms over the radial grid 
of log(p ; /p(r,)) with pj the density estimate at grid point rj 
and p the parametric fitting function; this is identical to how 
most observers define the residual. Below we state the rms 
deviation between the "measured" profile and the best-fitting 
parametric model in terms of magnitud es, denoted by A/u. 

We followed the practice in Navarro et ah (120041) of only 
constructing density estimates in the radial range r conv < r < 
7"2007 where r coml is the radius beyond which the halo mass 
distribution is considered robust to errors or approximations 
associated with the simulations (particle softening, relaxation 
etc.) and r2oo is the virial radius, i.e. the radius within which 
the mean densit y contrast is 2 00 times the critical density. Ta- 
ble 2 o flNavarro et alJ (120041) gives values of r2oo and r conv for 
all halo models. 

3. DARK MATTER HALOS AS SERSIC MODELS 

With few exceptions, modelling of the luminosity profiles 
of galaxies is done in projected space. We therefore began by 
analyzing the surface density profiles of the dark halos. Non- 
parametric estimates of T.(R) for the 19 halos are shown in 
Figure la, and Figure lb plots the deviations from the best- 
fitting Sersic model, equation Q; Table 1 gives the best- 
fitting n and A/j. The mean Sersic index is 3.11 ±0.49 (D), 
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3.00 ±0.17 (G), 2.38 ±0.24 (C), possibly indicating a (weak) 
trend toward decreasing curvature (lower n) in the profiles of 
halos of increasing mass. 

The Au va lues average 0.043 CP), 0.048 (G), 0.04 8 (C). For 
comparison. lCaon. Ca naccioli & D'Onofrio ( 1993) fin d Ay » 
0.05 in a sample of 45 E and SO galaxies, and Truiill oet al.l 
(2004) find a mean A/j of 0.09 in a sample of 12 elliptical 
galaxies without cores. The radial range over which luminous 
galaxies are fit varies from ~ 1.5 to ~ 3.5 decades, compara- 
ble on average with the ~ 2 decades characterizing our dark- 
matter halos. While the noise properties are different for the 
two types of data, the particle numbers in our halo models are 
small enough (~ 10 6 ) to contribute nonnegligibly to A/j. We 
conclude that the Sersic law fits dark halos as well as, and 
possibly even better than, it fits luminous galaxies. 

The residuals in Figure lb appear to show some structure. 
We will return in a future paper to the question of whether 
a modification of the Sersic law might reduce these residuals 
even further. 

4. WHICH FUNCTION FITS THE SPACE DENSITY BEST? 

Navarro et al. (2004) showed that equation (|4} provides a 
good fit to the spatial density profiles of dark halos. We 
showed above (§3) that the Sersic law Q is a good fit to the 
surface density profiles of dark halos. An obvious inference 
is that a deprojected Sersic law should provide a good fit to 
the space density. Here we ask which function - equation (0}, 
or a deprojected Sersic law - gives a better fit to p(r). We also 
consider the quality of fit of another three-parameter function, 
the generalized NFW profile presented in equation (|5}. 

When fitting deprojected Sersic profiles to the dark halos, 
we define n c i to be the Sersic index of the projected function; 
hence n e / should be close to the index n derived when fitting a 
Sersic law to the surface density (and the two would be equal 
if the halo's surface density were precisely described bySer- 
sic's law). When reporting fits to p(r) with equation (@J we 
define n = oc , with a the shape parameter of equation 0. 

The results are shown in Figure 2 and Table 1 . Mean val- 
ues of A/j for the three fitting functions (deprojected Sersic, 
eq. ®, generalized NFW) are (0.043,0.047,0.073) for the 
dwarf halos, (0.054,0.055,0.060) for the galaxy halos, and 
(0.061,0.063,0.060) for the cluster halos. Thus, the two Ser- 
sic functions are almost indistinguishable in terms of their 
goodness of fit: at least over the radial range available, a 
deprojected Sersic profile with index 2.5 <^ n c / <, 3.5 can be 
well approximated by a Sersic profile with n in the range 
4.5 ^ n <, 7.5. Both functions provide a significantly better 
fit to p(r) than the generalized NFW profile in the case of 
the dwarf halos, and the two Sersic functions perform at least 
slightly better than NFW for the galaxy halos. No single func- 
tion is preferred when fitting p(r) for the cluster halos. 1 

Another way to compare the halo density profiles with Ser- 
sic's law is via the radial dependence of the slope. Figure 
3 shows nonparametric estimates of the logarithmic slope, 
d log p jd log r, for the dark halos; slopes were computed via 
direct differentiation of the kernel density estimates, using a 
larger kernel width to compensate for the greater noise gen- 
erated by the differentiation. Equation (0 predicts a straight 

1 Also of interest are the mean y- values in the fits to the generalized NFW 
profile. We find (y) = (1.24, 1.26, 1.17) for dwarf, galaxy and cluster halos 
respectively. We note that these are significantly shallower than the steep 
inner slope (y = 1.5) proposed by Moore et al. (1999) and, as discussed 
by Navarro et al. (2004), are best interpreted as upper limits to the inner 
aymptotic behavior of the profile. 
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FIG. 2. — (a) Nonparametric estimates of the space density of the 19 dark 
halos. Vertical normalization is arbitrary, (b-d) Deviations in magnitudes of 
three parametric models from p(r): (b) deprojected Sersic model; (c) equa- 
tion |4}; (d) generalized NFW model, equation |5j. Best-fit parameters are 
given in Table 1 . 



line on this plot. That is a reasonable description of Figure 
3. The value of d log p/d log r in the G and C halos reaches 
~ — 1 at the innermost radii, consistent with the asymptotic 
power-law inner behavior of an NFW profile. No obvious 
convergence to a power law (constant logarithmic slope) is 
seen in Figure 3, and it is likely that simulations of improved 
resolultion may l ead to even shallower slopes at smaller radii, 
as pointed out by Navarro et a il d200l . 

5. WHAT DOES IT MEAN? 

Figure 4 shows Sersic's n (derived from fits to the surface 
density) as a function of mass for our dark halos and for a sam- 
ple of early-type galaxies. There is overlap at M « 1O 1O M0, 
the mass characteristic of "dwarf" halos and giant ellipticals. 
However the galaxies exihibit a much wider range of n values, 
extending to n < 0.5 in the case of dwarf ellipticals. A natural 
interpretation is that n is determined by the degree to which 
(dissipationless) merging has dominated the evolution. The 
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FIG. 3. — Nonparametric estimates of the logarithmic derivative of the 
space density for the 19 halo models. 
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FIG. 4. — Sersic index (derived from fits to the surface 
density) versus mass for gal axies (open circ l es) an d dark halos. 
Galaxy points are taken from [Bmggeli & Jerien 1 1998); StiaveU^^jJ 
l200ll): IGraham & Guzman] <20q3D:ICaon. Capaccioli & D'Onofrid 119931) * 
D'Onofrio. CaDaccioli & Caon 1 1994,). Halo masses are M200 from 
Navi"rro"etiuTl2'ob4l) . Galaxy masses were computed from total luminosi- 
ties assuming the Magorrian et al. (1998) mass-to-light ratio, with Ho = 70 
km s~ 1 Mpc~ 1 . 



nearly exponential (ra w 1) profiles of dE galaxies are similar 
to those of disk galaxies, suggesting that dissipation played a 
critical role in their formation. Luminous ellipticals are the 
end products of many mergers, the most recent of which are 
likely to have been gas-poor, and have de Vaucouleurs-like 
profiles (n w 4). This view is supported by numerical sim- 
ulations (Scannapieco & Tissera 2002; Eliche-Moral et al. 
2005) that show how exponential profiles are converted into 
de Vaucouleurs-like profiles via repeated mergers. 

A thornier question is: Why should a law like Sersic's 
fit dark or luminous spheroids in the first place? Ser- 
sic's law with 2 <, n <, 4 has an energy distribution that is 
roughly Boltzmann, N(E)dE ~ e$ E dE, and it is sometimes 
loosely argued that this "maximum-entropy" state is a re- 
sult of t he mixing that accompanies violent relaxation or 
merging jBinne\ll982tlMerritt. Tremaine & Johnstone! 19891 
ICiottilll991l) . With regard to dark halos, Taylor & Navarro 
(2001) have shown that the dependence of phase-space den- 
sity on radius is well approximated by a power law whose 
corresponding inner density profile has the shallowest slope. 
This can again be interpreted as an indication that the halos 
are well mixed. While our study does not shed a great deal of 
light on this question, it does suggest that the scale-free prop- 
erty of Sersic's law, (ilnp/iilnr °< r a , is the feature that links 
dark and luminous spheroids and that this property may be a 
hallmark of systems that form via gravitational clustering. 

We have shown that the fitting function that best describes 
luminous galaxies, the Sersic law, is an equally good fit to 
dark halos. We have not shown that the Sersic law is a good 
fit in an absolute sense to either sort of system. But given that 
dark and luminous density profiles are not pure power laws, a 
three-parameter law like Sersic's is as parsimonious a descrip- 
tion as one can reasonably expect. Future work should explore 
whether other, three-parameter fitting functions can describe 
dark and/or luminous systems better than Sersic's law. 

We thank Alister Graham for making available the galaxy 
data in Figure 4 and for informative discussions. DM was 
supported by grants AST-0206031, AST-0420920 and AST- 
04375 19 from the NSF and grant NNG04GJ48G from NASA. 
JFN acknowledges support from the Alexander von Humboldt 
Foundation. 
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